pvi_v1 <- lm(violence ~ days_scaled * treatment + pvi_state, stan_violence, weights = weight)
summary(pvi_v1)
confint(pvi_v1)

pvi_v2 <- lm(violence ~ days_scaled * treatment * pvi_state, stan_violence, weights = weight)
summary(pvi_v2)
car::Anova(pvi_v2)

pvi_n1 <- lm(norms ~ days_scaled * treatment + pvi_state, stan_norms, weights = weight)
summary(pvi_n1)
confint(pvi_n1)

pvi_n2 <- lm(norms ~ days_scaled * treatment * pvi_state, stan_norms, weights = weight)
summary(pvi_n2)
car::Anova(pvi_n2)

pvi_a1 <- lm(affpol ~ days_scaled * treatment + pvi_state, stan_affpol, weights = weight)
summary(pvi_a1)
confint(pvi_a1)

pvi_a2 <- lm(affpol ~ days_scaled * treatment * pvi_state, stan_affpol, weights = weight)
summary(pvi_a2)
car::Anova(pvi_a2)

stargazer::stargazer(pvi_v1, pvi_v2, pvi_n1, pvi_n2, pvi_a1, pvi_a2,
                     type = 'latex',
                     title="Results by PVI",
                     omit.stat = c("ser", "adj.rsq", 'rsq'), 
                     font.size = "scriptsize", column.sep.width = "6pt",
                     df = F, no.space = T,
                     dep.var.labels = c(
                       "Violence",
                       "Norms", 
                       "Aff Pol"))
